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Abstract 

A group of transformations changing the phases of the elements of 
the single-particle density matrix, but leaving unchanged the predic- 
tions for identical particles concerning the momentum distributions, 
momentum correlations etc., is identified. Its implications for the de- 
terminations of the interaction regions from studies of Bose-Einstein 
correlations are discussed. 

PACS numbers 25.75. Gz, 13.65.+i 
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1 Introduction 

Much work is being done on Bose-Einstein and Fermi-Dirac correlations in 
multiparticle production processes. The purpose in most cases is to deter- 
mine the features of the interaction region i.e. of the region where hadrons 
are produced. Let us suppose that function p(x) is the density of some kind of 
hadrons, say of negative pions, just after they have been produced. Knowing 
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this function one would know the geometry of the interaction region. Ac- 
cording to textbook formulae p(x) is proportional to the diagonal elements 
of the density matrix^ in the coordinate representation pc(x;x). In multiple 
particle production processes, however, these matrix elements cannot be di- 
rectly measured - there is no ruler to measure the interaction region. One 
measures momenta and uses the relation between the density matrices in the 
coordinate and the momentum representations 

Pc(x;x) = J ^^|^p(p,;P2)e-^ (1) 

where 

q = Pi - P2 (2) 

Thus, in order to determine the shape of the interaction region it is nec- 
essary to know both the diagonal and the out-of-diagonal elements of the 
density matrix in the momentum representation. The diagonal elements are 
easily obtained, because they are proportional to the density of the single 
particle momentum distribution ^^(p). The only thing one can do to get 
information about the out-of-diagonal elements is to measure the two-body, 
three-body etc. momentum distributions. In general it is not possible to 
express these distributions in terms of the single particle density matrix. 
Making, however, the usual assumption that the fc-particle density matrix 
can be approximated by the symmetrized product of single particle density 
matrices (see, e.g., one can determine from the momentum distri- 

butions of up to k particles the functions 

^ [p(Pi; P2)p(P2; Ps) • • • p(Pfe; Pi)] . (3) 

where 3? stand for real part of. Thus from single particle momentum distribu- 
tions one gets the diagonal elements p(pi; Pi) as already mentioned. Adding 
the information about two-body momentum distribution one finds further 
p(Pij P2)p(P2; Pi) = |p(Pij P2)P ^ the absolute values of the out-of-diagonal 
matrix elements, but not their phases. Measuring more-particle momentum 
distributions one finds some further information abut the phases, but not all 
of it. It turns out that it is possible to include in the single particle density 

^It is simplest to use the interaction picture. Then the density matrix after freeze-out 
is time-independent 
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matrix elements an additional phase factor without changing the predicted 
momentum distributions. The choice of this factor, which is unconstrained 
by the data, has a dramatic effect on the size and shape of the interaction re- 
gion which one deduces from the momentum measurements. It is the purpose 
of the present paper to investigate this effect in a systematic way. 

In the next section the phase ambiguity is explicitly formulated. Its effect 
on the size and shape of the interaction region is discussed in Section 3. Two 
specific examples are presented in Sections 4 and 5. Our conclusions are 
summarized in the last section. 

2 Group of transformations leaving the mo- 
mentum distributions invariant 

It is easily checked that making the replacement 

p(Pi; P2) - p^(Pi; P2) = e^[*(p^)-'^(p^)V(Pi; P2), (4) 

where 0(p) is an arbitrary, real- valued^ , function, one leaves all the observ- 
ables Q unchanged. Thus the measured momentum distributions are in- 
variant with respect to the transformation given by (0)). In other words, 
the measured momentum distributions yield no argument for or against any 
choice of 0(p). 

Before analyzing the effect of this phase ambiguity on the determination 
of the interaction region let us recall the standard procedure. 

In order to obtain information on the interaction region from the Bose- 
Einstein correlations in multiple particle production one usually applies the 
formula (cf. e.g. the review jl]) 

P(Pi;P2) = I d'XS{K,X)e"^'', (5) 

where 

K = ^{p,+P2), X = ^{XI+X2) (6) 

^When (f> is real-valued, matrix is hermitian and has trace one, thus it can be 
interpreted as a single particle density matrix. 
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and S{K, X) is known as the emission function. It follows from the her- 
miticity of the density matrix that the emission function is real- valued. It is 
constructed by analogy with the Wigner function and is interpreted as the 
position and momentum distribution of the produced hadrons As 
is well known, this interpretation cannot be exact^. We will accept it here, 
however, without further discussion. 

Another well-known difficulty (see e.g. [S]) in the studies of the interaction 
region is that the Fourier transform (jS)) cannot be inverted. The reason is 
that the density matrix on the left-hand side is known only for the on-mass- 
shell values of the four- momenta pi,P2- In general, there is an infinite variety 
of different emission functions which correspond to a given density matrix. 
In order to avoid getting involved with this difficulty we limit our discussion 
to models, where the emission function can be unambiguously expressed by 
the Wigner function. Two such cases have been much discussed. When all 
the identical particles are produced simultaneously at some time t = 0, one 
can choose 



where W is the Wigner function corresponding to the time-independent (in- 
teraction representation) density matrix in the momentum representation. 
When the distribution of longitudinal momenta is weakly correlated with the 
distribution of transverse momenta, one can consider the transverse emission 
function S't(Kt, as the Wigner function corresponding to the transverse 
factor Pt(Pti) PT2) in the density matrix We will thus use the relation 



where p^(K,q) = p^{Pi',P2), dq = (iq/(27r)' and dq = d'-q with / being the 
number of dimensions considered ( the vectors are two- or three-dimensional 
depending on the model). The Wigner function VF^(K, X), given by (jHI), will 
be used to define the interaction region. 

•^Taken literally, it contradicts the uncertainty principle. The solution is also well 
known: one has to perform an adequate "smearing" of S{K, X) bringing it into agreement 
with the requirements of quantum mechanics. 

^For an example see Section 5. 




(7) 
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3 The phase and the moments of the X dis- 
tribution 



In order to study the effects of the full group of transformations (j3)), let 
us first discuss the moments of the distribution of X. The i^'-dependent 
averages, for a given choice of the phase (j), are denoted by 

and the full averages by 

{{g))^ = j dKdXW^{K,X)giX). (10) 

Since in the latter average the Wigner function has been integrated over 
K, there is no problem with the uncertainty principle^. It is convenient to 
introduce the notation 



(X)^=o(K) =ro(K); {X')^=o = ^o(K); (11) 
((X))^=o =ro; ((X2))<^=o =Rl (12) 

Using the definition (jH)) and the identity (one dimension, n = 0, 1, . . .) 

dxx^'e'"'' = {-^T-Q^ j dxe''^'' = 27r(-z)"— 5(g). (13) 

one finds that the denominator on the right-hand side of ^ equals P(f,(K, 0), 
and 

'Tl / i^TT. \ 

which is easily generalized to 

'Tl I i^'M \ 

<n^rwK)-^fip^p*(K,q)) , (15) 



q=0 



where n = uj . 



^There may be, however, problems with the interpretation of the averages (g(X))(K). 
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When evaluating these averages we will need derivatives of the phase 
factor 



These formulae allow to calculate all moments ((X")) in terms of (the 
integrals of) the derivatives of the single particle density matrix. This gives 
the full information about the size and shape of the interaction region. 

4 Position-momentum correlations and the size 
of the interaction region 

We will now discuss in detail two parameters characterizing the interaction 
region: (i) the average position ((X)) and (ii) the averaged square of the size 



[V (0(Pi) - 0(P2))]q=o = V0(K); [V^ (0(pi) - 0(P2))]^^o = (17) 



0(Pi) - 0(P2) = 0(K + ^) - 0(K - ^) = 2 sinh(^)0(K), (16) 




one finds 



(X)^(K) 
(X^)<^(K) 



ro(K)- V</.(K), 

RliK) - 2ro(K) ■ V0(K) + (V0(K))' 



(18) 
(19) 



and 



((X)), 
((X^))<^ 



ro - (V0(K)) 

i?2-2(ro(K)-V0(K)) + ((V0(K))^), 



(20) 
(21) 



where the notation 




(22) 



^They are special cases of 
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has been used. 

Function (X)^(K) can be interpreted as a measure of the position-momentum 
correlations. According to formula (jU the phase is fixed, when the function 
(j) is given. On the other hand, as seen from formula (fTHjl . 0(p) can be also 
fixed, up to an uninteresting constant, by specifying the position momentum 
correlation (X)(K). 

The variance of X at fixed K is of particular interest because it gives the 
HBT radius, as determined in the usual way from the two-particle distribution 
in the Gaussian approximation. This can be seen as follows. First, we observe 
that it is invariant, i.e. does not depend on the choice of 0(p) : 

(X2)^(K) - ((X)^(K))^ = RliK) - rl{K). (23) 

It can be therefore calculated independently of the choice of 0(p). On the 
other hand, in the Gaussian approximation 

p(K, q) = p(K, q = 0)e-'i'RHBTi^)/^ (24) 

where Rf^^j-lK.) is defined by this relation. Formula (|14j) implies for this 
approximation (X)(K) = 0, {Xf)(K) = i?^^y(K). According to (j2Sl) we 
thus have 

dRlM^) = (X^)(K) - ((X)(K))^ = RliK) - rl{K). (25) 

where c? = 2, 3 is the number of dimensions. 

Note, however, that the variance of X at fixed K (and thus also the 
measured HBT radius) is not the physical size of the interaction region at 
fixed K. It is only an auxiliary concept. Indeed, according to Heisenberg's 
uncertainty principle fixing K implies that the true variance of X becomes 
infinite (cf. footnote 3). 

The true size of the interaction region, as derived from the Wigner func- 
tion, is instead determined by the full average, given by 



i?Le = ((X^))0-((X))| = i?^-r^ + 
+ ([V0(K)]2) - ((V0(K)))2 - 2((ro(K)V0(K))K - ro(V0(K))) (26) 

One sees that this quantity explicitly depends on </)(p) and so does the in- 
terpretation of the HBT measurements. 
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The three pairs of terms on the right-hand side are easily identified as: 
the variance at fixed = 0, the variance of V0, and twice the unnormahzed 
correlation coefficient for ro(K) and V(/>. The variance of increases the 
observed variance. The correlation may have either sign. In particular, it 
may be negative and so large in absolute value that it reduces the observed 
variance below the variance at = 0. This is not surprising. Suppose that 
the variance for some density matrix po is smaller than for p^. Then one 
could start with and go over to the original po by introducing a suitable 
phase. This change of phase reduces the variance. 

Equating to zero the (variational) derivative of the right-hand side of 
equation (j26j) with respect to V0 one finds the condition for the variance of 
X (i.e. the size of the interaction region) to be minimal. It turns out that 
the minimum is achieved when is chosen such that the average of X at 
fixed K does not depend on K: 

(X)^(K) = ro(K) - V0(K) = a, (27) 

where a is a constant vector. In other words, the minimal size is obtained 
when the phase is chosen such that there is no position-momentum corre- 
lation. It is clear from ()27|) that one can always find which satisfies this 
condition. 

Substituting (|77j) into (fTH|) one finds that the minimal size of the system 
(for the given momentum distribution) is 

[((X^)),-((X))J]_ = Rl-{vl{K))=<Rl^^{K)>. (28) 

For a generic 0(p) the interaction regions at various K's are shifted with re- 
spect to each other. This increases the overall, integrated over K, interaction 
region. Condition (X)0(K) = a puts all these partial integration regions, as 
well as possible, on top of each other. This minimizes the overall, calculated 
interaction region. The minimum is just the average over K of expression 
(OHj) . We thus conclude that the HBT radius, averaged over K, corresponds 
to the minimal possible size of the interaction region. 

In order to find the moments up to order 2m it is enough to know the 
odd-order derivative of up to V^"'~^0. Thus (X) and (XiXj) can be calcu- 
lated when V0 is known. Knowing moreover V^0 one can calculate all the 
moments up to fourth order. In general it is not possible to calculate the 
Wigner function PV^ in closed form. We will discuss, however, two instructive 
special cases. 
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5 Example I 

We shall consider only the transverse variables. Choosing 

0(p) = a + Y, b^Pi + Yl ^^Pi ' (29) 

i i 

with a,bi,Ci constant {i = x,y), one obtains 

0(Pl) - 0(P2) = Yl l^^^i + = Y (30) 

i i 

According to formula (|18p this is the most general case where the phase 
4>{pi) — <f){p2) is linear in q. In this case, from (0), we can write 

P^(Pi;P2) = / rf'X5(K,X)e^''(^+^). (31) 

Changing the integration variables X into X + V one finds the emission 
function corresponding to the density matrix 

^^(K,X) = S(K,X-V)) (32) 
Thus, the interaction region is shifted in space-time by V: 

X,-^X,-bi- 2c,K,. (33) 
This formula shows that introducing the phase ()29|) may have two effects: 

• A rigid shift of the interaction region. 

• A correlation between the momentum and the space-time position of 
the emission point of the particle. 

The first point is well-known and not really disturbing. It has been discussed 
by many authors. The second has some relation to the notion of the ho- 
mogeneity region jUl |7j , but to the best of our knowledge its relation to the 
phase of the density matrix has not been pointed out. 

In order to show that these considerations are not purely academic we 
will now consider a specific model [3 UHl HI] where the emission function is 
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derived from a physical picture. The main assumptions, at the classical level, 
are 



for 



/i = 0, 1,2,3, 



(34) 



X 



T 



(35) 



where r is a constant. Assumptions (j34p and (j35|) for the temporal and 
longitudinal components were introduced by Bjorken [^1 and Gottfried ^3] 
in their boost invariant description of the multiple particle production pro- 
cesses. The full assumption (jSH) was introduced by Csorgo and Zimanyi 
who used it to explain why in e'^e~ annihilations the correlation function 
depends on only. The motivation in ^ was to explain why the estimated 
radii of the interaction region decrease rapidly with increasing particle mass. 
Let us add the definition of the transverse mass 



Note that, while for a particle of given mass m; 



(36) 

+ pf. depends on 

the transverse momenta, here is not fixed and thus M|, depends on the 
longitudinal and temporal components of momentum. Assumptions ()34|) and 
fjH^ imply 



A 



Mt 



T 



(37) 



It is easy to see qualitatively, why in this model the measured radius of the 
interaction region is a decreasing function of the particle mass. Let us make 
the crude assumption that the momentum distributions are the same for all 
the particles. Then the distribution of MtX is universal and since M-r is 
bigger for heavy particles, the corresponding distribution in X is narrower. 
Let us consider the transverse part of the emission function proposed in 

uni: 



T 



exp 



'24 



(38) 



It can be interpreted as a Wigner function^ and thus our previous discussion 



^At fixed Mt 
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is applicable to it. This Wigner function can be rewritten in the form 



St = exp 



(X 



T 



2RI 



where 



2m 



(39) 



R(j, 
Rd 



l + fi^ St 



Tt 



with 



liR^ 



TtMj 

t6t 



(40) 

(41) 
(42) 



(43) 



St and (pT are of the form discussed in the preceding section^. We thus 
conclude that the measured HBT radius is R^, whereas the true size of the 
interaction region, according to (jSHI), is 



true 



(44) 



Thus, while R^ decreases with increasing transverse mass Mt, Rtme does 
not depend on the transverse mass. This explains why the measured HBT 
radii may be very small at large particle (transverse) mass, while the actual 
emission region is mass independent. 



6 Example II 

As an instructive, more complicated, example of the phase we now consider 

i 

^Because of the presence of Mt this is an approximation, vahd when only two-body 
symmetrization is used and/or when the variance of Mt can be neglected. 
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where i = x,y (we consider only transverse dimensions). In order to keep 
the discussion simple we will start with the Gaussian density matrix 



P(Pl,P2) 



27rA2 



With this choice we obtain 



(46) 



(47) 



where we have used the identity pi + pip2 + pi = S-ft'^ + ^q^ valid for both 
the X and the y components. Using the formulae of the previous section we 
find from ()46|) and ()47p the average shift of the X distribution: 



< X, > {K, 



(48) 



where j = x,y. 

From relation (|47p we obtain 



SJK,X) 



47r2' 



,iqX 



Pi,P2) = s^^{K:,,x)sy^{Ky,y), 



(49) 



with 



^ 27r 



exp 



(50) 

The integral can be evaluated in terms of the Airy function. To simplify 
notation, from now on we drop the subscript j. 
Introducing a new variable z by 



where 



we obtain 



q = a(z uj 

^ ^ 2 



UJ = aR 



(51) 
(52) 



ae 2a2 



AV27r 



exp 



-z^ +iAz + B 



dz 
2^' 



(53) 



where 
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A = a[X- <X > (K)] + LO^A = aX + 0:^4, (54) 



B = ^ 
2 



2 



aX + ujyQ (55) 



The integral can be expressed in terms of the Airy function (cf. [15^ for 
an equivalent formula) 

1 r+oo 

Ai{u) = / ei*'+'"*rft (56) 

The result is 

s^{K,X) = ^^e^e^M{A) (57) 

One sees from this formula that, for fixed K, the variation of X, i.e. 
spread of X around the average {X){K) is determined by two parameters: a 
and u. From the known asymptotic expansion of Ai{z) one can deduce that 
in the limit a ^ oo one recovers the original Gaussian distribution. In the 
limit a ^ the distribution becomes singular. Formally one obtains 



sJK,X) —Ai(aX). (58) 

and the averages of X and do not exist. 

The Airy function is negative for some values of A, showing explicitly 
that the function s^{K^ X) cannot be literally interpreted as a distribution 
function. This is illustrated in Fig. 1, where s^{K,X) is plotted versus aX 
for K = Q and = 0, 1, 2. For = it is an Airy function. With increasing 
uj the maximum moves towards X = 0, the oscillations get dumped and 
already at = 2 the curve is almost a Gaussian. 

The formulae from Section 3 yield 

{xj) = ^; a^{xj) = R^ + 2{xj)^ (59) 

The net result is a shift of the interaction region and an increase of its size. 
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Figure 1: Emission function s^{K,X) versus aX. Dotted line: = 0, dashed 
line: a; = 1, full line: lo = 2. 

7 Conclusions 

For models where particle production is uncorrelated except for the Bose- 
Einstein correlations, we have identified a group of transformations of the 
phase of the single particle density matrix which leave all the the single- and 
multiparticle momentum distributions invariant. The effect of the resulting 
uncertainty of the size and shape of the interaction region as deduced from 
the measured momentum distributions is systematically discussed. Explicit 
formulae taking into account this uncertainty are derived for the moments of 
the X distribution. It is shown that the phase ambiguity can have a dramatic 
effect on the parameters characterizing the interaction region. Two explicit 
examples are discussed in detail. 
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